Inhomogeneous ground state and the coexistence of two length scales near phase 

transitions in real solids 
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Real crystals almost unavoidably contain a finite density 
of dislocations. We show that this generic type of long-range 
correlated disorder leads to a breakdown of the conventional 
scenario of critical behavior and standard renormalization 
group techniques based on the existence of a simple, homoge- 
neous ground state. This breakdown is due to the appearance 
of an inhomogeneous ground state that changes the character 
of the phase transition to that of a percolative phenomenon. 
This scenario leads to a natural explanation for the appear- 
ance of two length scales in recent high resolution small-angle 
scattering experiments near magnetic and structural phase 
transitions. 
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The modern theory of phase transitions (PTs) in pure 
and weakly disordered crystals is based on the scaling 
hypothesis and the renormalization group (RG). It has 
been confirmed by a large number of experiments and 
numerical simulations. 

Therefore, the results of some recent high resolution 
X-ray and neutron scattering experiments for various 
crystals near structural and magnetic PTs were quite 
unexpected 0] : in contrast to the predictions of conven- 
tional scaling theory two distinct large scales were ob- 
served. The temperature dependence of the smaller scale 
was found to be broadly consistent with the results of 
the conventional theory for the correlation scale of ther- 
mal fluctuations. Quite recently explicit experimental 
evidence for Verneuil grown samples of SrTiOs Q con- 
firmed earlier speculations that the second, larger scale 
is connected with the presence of defects in a very large 
region near the surface. 

All experiments that found two length scales used crys- 
tals which were considered to be of good quality. Nev- 
ertheless, surface preparation creates a significant den- 
sity of defects in a surface layer with a thickness up to 
several /x [0-0. Based on this observation, Altarelli ct 
al. suggested [S] that most of the defects are actually dis- 
location dipoles which induce elastic strains with a pair 
correlator G{r) ~ |r|~^- They tried to explain the criti- 
cal behavior (CB) in the disordered layer as that of the 
Weinrib-Halperin (WH) model The coexistence of 
two scales was then interpreted as a simple superposi- 



tion effect, with scattering from the larger scale in the 
disordered layer (where the CB is governed by the WH 
fixed point) and from the smaller scale (corresponding to 
conventional CB) in the bulk of the sample. However, 
later experiments using very thin Holmium films found 
clear evidence that near the magnetic PT temperature 
both scales coexist in the same volume fraction of the 
sample 0. 

By construction, the conventional RG procedure is not 
able to explain the simultaneous existence of two length 
scales. Not only its detailed value and temperature de- 
pendence, but even the principal theoretical basis for the 
origin of the second length scale appears to be unclear. 

In this Letter we point out the shortcomings of the RG 
procedure for the treatment of systems with weak, but 
long-range (LR) correlated disorder, the generic situa- 
tion in crystals with topological defects like dislocations. 
We show that the global scenario of the PT differs from 
that of models with short-range correlated disorder: in 
addition to the thermal correlation length (which does 
not diverge, but saturates near T^) even for very small 
dislocation densities a second length scale inevitably ap- 
pears. Its divergence determines the critical point. In 
the language of the RG, this scale is connected with 
the existence of a temperature dependent inhomogeneous 
ground state (GS). This GS is formed by stable nuclei of 
the ordered phase that appear even above the global Tc- 
We compare our results with experiments for some per- 
ovskites and rare earth magnets (REM) [|l|. 

We begin with a brief review of the standard RG treat- 
ment of the CB of weakly disordered crystals. The Hamil- 
tonian of the simplest model of this class, for an n- 
component order parameter ipi, i = I ... n, is given by 



d'r^AT{r)^Ur) (1) 



with the Ginzburg-Landau Hamiltonian Hgl for a pure 
crystal. 

The disorder variable AT(r) is usually assumed to be 
Gaussian distributed, so it is completely characterized by 
the pair correlator 



G(r) = (AT(O)AT(r)). 



(2) 



This can be described by two parameters, a defect 
strength rj = G(0) and a single, finite correlation length 
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Rq, so G(r) — > for |r| :§> Rq. For actual RG calcula- 
tions a further simplification is usually to ignore the scale 
i?o and set 

Gir)^uoS{r); uq = R'^. (3) 

This is selfconsistently justified with the argument that 
close to Tc any finite scales like i?o are irrelevant in com- 
parison to the diverging scale of critical fiuctuations, 

However, the RG procedure is a description of the criti- 
cal fluctuations around the GS of a system, and its results 
are valid only if the GS is determined correctly. Usually, 
this GS is assumed to be homogeneous; but in fact it 
has to be determined using either a Landau type saddle 
point equation for a particular disorder configuration (see 
below), or the saddle point equation for the disorder av- 
eraged free energy. For uncorrelated Gaussian disorder, 
the latter can be written as 

(r-A)^f (r) = K-Ao) (^^(r))', (4) 

with the pure vertex Aq. 

If Aq ^ uq this equation only has the trivial solu- 
tion ipf^\r) = for r > 0. But for stronger disorder. 
Mo > Mc = ^Ao with ^ w 1 nontrivial solutions appear 
, and we have to take into account the presence of a 
nontrivial GS. The standard RG procedure may be used 
only if Mo ^ Ao to describe "weak" disorder in the con- 
ventional RG sense In the opposite case mo > Uc 
a RG calculation has to account for the presence of an 
inhomogeneous GS. 

This leads to the question what kind of GS exists 
in "high quality" samples of real solids. The disorder 
strength is small, so a change of the CB can be ex- 
pected only when the correlation scale Rq is so large that 
the approximation of (H) by is no longer valid. 

Several types of defects lead to LR correlated random 
temperature-like disorder. Some special cases, like in- 
homogeneously distributed point defects have been dis- 
cussed before But there is a much more generic 
class of such defects: topological defects like dislocations 
in crystals. Models with dislocations require a more de- 
tailed analysis of many different situations, with varying 
types of dislocation ensembles and transition scenarios. 

We begin our analysis of such crystals with dislocations 
based on the Landau expansion of the free energy density 
/, 

/ = Hgl W,} + q.jki 4? (r) ViVj- (5) 

The tensor qijki describes the coupling between the 
order parameter and the elastic strain introduced by the 
quenched disorder, which is given by the sum over con- 
tributions from all dislocations p^ , 

a,'i 



where 4>i is the angle between the component of r — 
perpendicular to the dislocation and the Burgers vector 
ha of a dislocation of type a at Tj, and the sum runs over 
all dislocations. 

To simplify the equations, we will limit the discussion 
to the case of a scalar order parameter ip and an isotropic 
elastic medium |l3|. The GS for a proper RG treatment 
can be identified minimizing the free energy (|^). This 
leads to the nonlinear stochastic equation 

-g Aip + a (r + I (r)) if + \ + D ip"^ = 

(7) 

The simplified coupling is q/a ~ (K/Tq) (dTc/dp) with 
the elastic modulus K , the Curie temperature To of the 
undisturbed system and pressure p. The inclusion of a 
tf^-term allows us to discuss the weakly first order struc- 
tural FT in perovskites. 

The statistical properties of the solutions of (Q) are 
determined by the random elastic strain e = e^^^\ Obvi- 
ously the coefficient of the linear term in (M) can become 
negative and allow a nonzero solution of ^), either due 
to the large contribution to e(r) of a single dislocation 
close to r, or because the contributions of several dislo- 
cations happen to add up to a large (negative) value of e 
at r. 

For randomly distributed dislocations loops with a ra- 
dius of curvature L much larger than the average dis- 
tance rd between loops the probability distribu- 
tion function for e in (^ is a Levy distribution [^ : 
Vie) = Vl^ie), in particular r{e) ~ e'^ for e > 1. 
The variance of this distribution is given by 

var e « (e(r,))' ln(L/r,) = (^^^ (8) 

with e{rci) ~ b/2TTrd and the corresponding temperature 
scale Ti. 

The pair correlation function for this distribution for 
r < L is G(r) = (e(O)e(r)) - na ln(L/r) with the dis- 
location density Ud — r^'^- It is determined by the con- 
tribution of many dislocations on scales rmin < r < L, 
while for distances r < rmin = i^'niL/rd))~'^^'^ the ran- 
dom value of e(r) is dominated by the contribution of the 
single nearest dislocation. 

As the elastic strains are strongest in the immediate 
vicinity of the dislocations, nonzero solutions of (^ may 
first appear there as the temperature is lowered from high 
values. In a first approximation, we can limit (|^) to in- 
clude only the strain induced by the closest dislocation i 
near r: 

~gA^ + a (r+^j^^^) ^ + X + D = 0. 
V a |r_L -Tiiy 

(9) 
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The temperature T„ at which nonzero solutions of 
this nonlinear, deterministic equation first appear can 
be found rigorously by an analysis of the linear part 
of this equation. The main results of this analysis are as 
follows. The reduced temperature t„ = (T„ — ro)/To and 
the characteristic scale of ordered nuclei i?(r„) are, 

T„ « — (g&)' + Ar,„ (10) 
ag 

R{Tn)^\^.f^^r,{rn). (11) 

Here Ar/i ^ r„ describes the temperature range of the 
hysteresis in a crystal without dislocations and a weakly 
first order PT like the perovskites that display the two 
length scales The second expression for i?(r„) re- 

flects the fact that nonzero solutions appear when the 
energy gain from the attractive potential balances the 
"kinetic" energy due to the gradient term in the free en- 
ergy. 

The above reasoning depends on two assumptions: the 
possibility to include only the simplest gradient term in 
(0) and (0), and the validity of the single defect approx- 
imation (|9|). 

Obviously, the former requirement holds if i?(T„) ^ a, 
where a is the lattice constant. The coefficient a can be 
written as a = Tq/C with the material's Curie constant 
C. In perovskites, for displacive PTs [|l|l C « 10"* ~ 
lO^if, so a = 0(10-2 - 10-3) and, indeed, i?(T„) « 
a/ a 3> a. The scattering experiments for SrTiOa also 
show that in the relevant temperature range ^ a, 
so the use of (0) and is justified. 

In REM, C < To Q, so a is large and (0) and (|) 
cannot be used at all. However, in both cases below r„ 
the radius R{t) grows much larger than Tc, so one can 
use the estimation 

R{r) « ^ {dTjdp) (12) 

in the temperature range where R{t) 3> ''c (''')• Obvi- 
ously the "divergence" of the scale R{t) is described by 
an exponent v = 1. The prefactor in ( [I4 ) has very simi- 
lar numerical values in perovskites and REM, which leads 
to scales of a few 10^ A in both classes In the tem- 
perature range where ( [T2| ) is valid, the integrated inten- 
sity of scattering from the nuclei of the ordered phase is 
I{T)^ndR^{T)^T-\ 

As the temperature is reduced the size and the num- 
ber of nuclei grow. To determine the global charac- 
ter of the PT one must then investigate the tempera- 
ture dependence of the disorder averaged GS correlator 
.g(r) = ((/3('''(0) ip'^^\r)) for r ^ 00. For this purpose it 
is crucial to analyze the state of "bridges" with t{v) < 
which link two nuclei that have appeared at some higher 
temperature [ p2[ . 

It is easy to show |^,^ that if on such a bridge r(r) < 
—T2 = g/ar"^^^ the correlation radius of local thermal 



fluctations is rc(r) ^ Tmin- In such bridges the energy 
of a domain wall (DW) is much larger than Tc, so no 
DWs will appear. If — T2 < r(r) < then rdr) > Tmin 
and DWs are easy to form. The probability for this is 
Prob{rc{r) > rmzn) « (t2/ti)^/2 = z. The sign of the 
order parameter in the nuclei connected by a bridge is the 
same with a probability 1 — 2 « 1 for z <^ 1. Of course, 
this sign will fluctuate in time as long as overlapping 
nuclei form finite clusters. 

As the temperature is decreased, these clusters will 
grow, and at Tc an infinite percolating cluster will ap- 
pear. On this cluster the sign of the order parameter can- 
not change anymore — this is the spontaneous symmetry 
breaking that defines the phase transition in the system. 
Therefore the global PT scenario is a percolative one. 
The condition 2: <C 1 can be rewritten as i?(T„) <C rmin 
using (ini). This turns out to be equivalent to the condi- 
tion under which the first nucleation of the ordered phase 
near dislocations can be described in the single defect ap- 
proximation (^. 

The smaller length scale seen in experiment corre- 
sponds to the fluctuations (f{r) = ip{r) — (y9(°)(r) around 
the inhomogeneous, temperature dependent GS i-p^^^ (r) 
determined by (|^) [Q. Assuming mean field exponents 
for simplicity, the width of the corresponding "broad 
component" in momentum space can be estimated as 

(Ag)2 = (r-2> = r' V{,, = f (r' - r)) dr' 

+ A_/!^|T'|7'(e=f(r'-r))dr' (13) 

with T = (T — Tq) /Tq. This implies that the line width 
remains finite at all temperatures, with a minimal value 
Agmin ~ '''c^i'^i)- This is reached (in the mean field ap- 
proximation, = 2) at p = 1/3. Because in three 
dimensions the percolation concentration pc < 1/3, this 
minimum lies slightly below the temperature Tc where 
the length scale connected with the narrow component of 
the scattering diverges, as observed in experiments 

If i?(r„) > rmin, neither the single defect approxima- 
tion nor the simple percolative picture of the PT hold, 
as the inhomogeneous GS (^'^^^(r) and thermal fluctua- 
tions both play an equally important role. In this case a 
completely new generalized RG scheme is needed. 

To summarize, we have analyzed the CB of solids with 
LR correlated quenched disorder taking into account the 
appearance of an inhomogeneous, temperature depen- 
dent GS. In real crystals the appearance of such a GS 
is a generic phenomenon due to the almost unevitable 
presence of dislocations. The appearance of this GS is 
important even for a qualitative discussion; in particular 
it leads to the breakdown of the conventional RG proce- 
dure. 

We have shown that if the condition i?(T„) ^ rmin 
holds the global PT scenario consists of two stages: at 
T = Tn nuclei of the low-temperature ordered phase ap- 
pear next to individual dislocation lines. At lower r 
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the nuclei grow and begin to overlap; the typical size 
of magnetically ordered domains is described by perco- 
lation theory. It is especially remarkable that this PT 
scenario becomes more likely when is reduced, be- 
cause rmin oc So generally speaking the two length 
scale phenomenon should be observable in scattering ex- 
periments in many crystals. Of course a high resolution 
setup with sufRcient sensibility as used in refs. is 
necessary, especially for 6etter samples with small n^, be- 
cause the volume occupied by nuclei of the ordered phase 

1/2 

becomes noticable at temperatures of order ti ^ nj 
(using and ([l^ ) one can write n = (udb'^y^^) for 
most structural PTs and the PT in REM. 

The universality of the percolation PT scenario due 
to LR correlated disorder can be illustrated by the the 
recent observation ||2^ of two scales in the neutron scat- 
tering data for the Invar alloy Fei_a;Nia; which is near a 
morphology boundary for x ~ 0.3. Contrary to earlier 
speculation ||] the description of the CB of real crystals 
very close to Tc must account for LR correlated disorder 
in addition to the uncorrelated disorder covered by the 
conventional theory. In fact, the asymptotic CB of "high 
quality" samples of perovskites and REM is the same as 
that of a "dirty" alloy with LR chemical disorder . 

This work was supported by SFB 237 (Disorder and 
Large Fluctuations) and by RFFI grants N16893 and 
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